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Near horizon geometries have been widely studied, and have found many applications. Certain 
static, near horizon geometries are now understood to be bulk duals to CFTs with static scale- 
invariant sources under the AdS/CFT correspondence. However, static near-horizon geometries 
aren’t just scale-invariant, they have extra ‘enhanced’ symmetry. This means that they can only be 
the bulk duals for a special class of static, scale-invariant sources that share this enhanced symmetry. 
The purpose of this paper is to consider bulk duals for more generic static, scale-invariant sources, 
without this extra symmetry. These solutions are quite different to near-horizon geometries. In 
place of the extremal horizon they have a null singularity. We find specific examples of such bulk 
geometries numerically for the cases of pure gravity, and for an abelian gauge field. 


I. INTRODUCTION 

The AdS/CFT correspondence [iHS] provides the opportunity to explore a particular class of strongly coupled field 
theories both in flat space, and on curved backgrounds. The dynamics of the CFTs that can be described by this 
correspondence is found, in the limit of large coupling and large numbers of degrees of freedom, by studying classical 
supergravity for spaces that are asymptotically locally AdSn+i x X, where X is some compact internal space. The 
conformal field theory lives on the n dimensional conformal boundary of this space. For given fixed choices of internal 
space X, we can always truncate to the Einstein’s equations, = Ag^"^ on an n -I- 1 dimensional asymptotically 
AdS subspace. We will take a bottom-up view, and also consider the addition of an abelian gauge-field, without 
worrying about whether it is a consistent truncation of supergravity, as is often done in AdS/CMT [4]. When we do 
this we use the Einstein-Maxwell equations 



= 0 . 

Restricting ourselves to these equations corresponds to examining a sector of dynamics in the large N, strong 
coupling limit, where the only vevs are the stress tensor, and a U{1) charge current which is sourced by the boundary 
value of the gauge field. This correspondence has been used in many contexts to probe strongly coupled field theories, 
both on flat and curved backgrounds. For a review of the use of AdS/CFT to study strongly coupled CFTs on curved 
backgrounds, see [5]. 

A particular subclass of the solutions to Q are the near-horizon geometries of extremal horizons. These are reviewed 
and classified in [S] . An important property of these solutions is their scaling symmetry, so, if they are asymptotically 
AdS, they can describe scale invariant states in a holographic CFT. Let us restrict to the static case, so that from the 
point of view of the CFT these describe static, scale-invariant states. As proven in [7], static near horizon geometries 
have an enhanced SO(2,l) isometry group^. This means they don’t describe the most generic CFT states in this 
symmetry class. 

For the remainder of this paper we specialize to 3-1-1 dimensions, but we expect that our results should generalize. 
In this case, all the static near-horizon geometries that solve 0 are known analytically [5]. The asymptotically AdS 


* Electronic address: a.hicklingl2@imperial.ac.uk 

^ Actually, what they show is that the near horizon geometry is a warped product of a maximally symmetric 2-d space, so the isometry 
group is the isometry group of either Minkowski, AdS, or dS. For the cases we’re interested in here, it’s AdS, with isometry group 
SO(2,l). 
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solutions have boundary metrics and electric potentials of the form^ 


T12 T 2 2 7 ; 2 2 1 i2 

g = —dt dv (X d(l) -^--t- ol dcj) 

, dt 
A = e —. 
r 


( 2 ) 


In 3 + 1 dimensions, we see that the isometry group is in fact S'0(2,1) x C/(l), with an additional C/(l) rotational 
symmetry in the (j) coordinate. 

We might, instead, want to consider the CFT with sources that explicitly break this SO{2, 1) x C/(l), either by 
introducing a boundary metric of the form 


a 


-dt^ + r ,\dr ,, 2 i ,2 

-- 'r2x{4>)—d(i) + a^d(lP, 


or by introducing an electric source with angular dependence 


A = F(</>) 


dt 


( 3 ) 

( 4 ) 


Introducing either of these sources breaks both the rotational symmetry, and part of the SO{2, 1), so that the only 
symmetry we are left with is the static, scale-invariant part. Thus the bulk dual clearly cannot be a near horizon 
geometry. The purpose of this paper is to consider these more generic bulks associated with static, scale-invariant 
states. As we will see, these solutions are quite different to near-horizon geometries, in particular the extremal horizon 
is replaced by a singular null surface. 

Finding these bulk duals involves finding 3-1-1 dimensional static scale-invariant solutions to 0 with negative 
cosmological constant. Once you impose the symmetries, this becomes an effectively two-dimensional problem that 
can be solved numerically, using the Harmonic Einstein method[8l [9]. We are able to solve for a range of choices of 
x(4') 8-nd V{(j)) on the boundary. For the case of pure gravity (no gauge field), these solutions were conjectured to 
exist in m- They were also constructed to linear order as a perturbation to AdS. Here we demonstrate that the full, 
non-linearized solutions exist by constructing explicit examples. 

The applicability of these results is not limited to scale-invariant boundary states. In [10] it was found that the near 
horizon geometries of extremal horizons could describe the large scale structure of holographic CFTs on geometries 
that are only asymptotically scale-invariant^. Specifically, it was found that if a static bulk contained an extremal 
horizon that met the boundary at infinity, then the near horizon geometry of this extremal horizon was determined by 
the geometry of null infinity on the boundary. In particular, static perturbations to the boundary geometry that fall 
off fast enough towards infinity would not affect the extremal horizon in the bulk, and so vacuum states would have 
the same large scale structure determined by this near-horizon geometry. These arguments were specifically for the 
case of pure gravity, however in m solutions have been found numerically where an electric potential is introduced 
on the boundary which is not scale-invariant, but is asymptotically of the form ([^. In this case the solutions found 
again had the near-horizon geometry which was determined solely by this asymptotic form. 

We would suggest that the appropriate way to describe the large scale structure of static CFT states on geometries 
that are asymptotically of the generalized form (|^, or where we have introduced electric potentials that are asymp¬ 
totically of the form 0. is through the bulks we describe here. In particular, these geometries would be suitable for 
the infra-red of the bulk duals to states on boundaries of this form. 


^ Note that the metrics in § are only determined up to a conformal factor. We can write g in the conformal frame 

g = —dt^ + dr^ + a^r^d(l)^, 

where we can interpret them as cones, or as fiat space when a = 1. The electric source is invariant under this change of frame. 
^ The conformal freedom means that asymptotically scale-invariant metrics include asymptotically locally flat metrics. 
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II. STATIC SCALE INVARIANT GEOMETRIES 

In Appendix we argue that a general static, scale-invariant, 2-1-1 dimensional boundary metric and electric 
potential can be written in the form 


rjf'^ fir 

9 = -^ + ^+ 

A = Vic^)-, 
r 

where the scaling symmetry is given by 


( 5 ) 


t —y ct 
r —>• er. 


( 6 ) 


Writing the metric in this form involves making a choice of coordinate frame. We can think of these static, scale 
invariant geometries as describing the large scale limit of boundary geometries with an asymptotic region of the form 

g = —dt^ + dr"^ -\- a^r'^dcj)^ — 2rax(0)d(^dr, (7) 

by Weyl scaling the metric by a factor of . In this frame 0. the surface r, t —>■ oo is null infinity, . If x((/<) = 0, 
it is an extremal horizon in <§ corresponding to the Poincare horizon of the AdS 2 factor. However, if xi'f) 7^ 0i this 
surface is singular. Pragmatically, this can be seen because the metric in (l7l cannot by transformed into Gaussian 
Null Coordinates, and this should be possible for any smooth hypersurface|l2|. This is the same singularity as we 
will see in the bulk of our solutions. Even if x{4’) = 0; turning on a non-constant V{(j)) will generate the singularity 
in the bulk. The nature of this singularity, and the resulting singularity in the bulk solutions, will be explained more 
physically in Section [V| 


A. Extremal Near Horizon Solutions 


On the other hand, when we consider boundaries with x{4') = 0) ^md constant V{(p) = e, there are non-singular 
bulk solutions that preserve the Symmetry <©■ The boundary geometry in these cases is a cone with opening angle 
determined by a. Alternatively, we can reinterpret this in the conformally compactified frame given in ([^. From this 
point of view it is AdS 2 x S^, and the parameter a now determines the relative sizes of the AdS 2 and Si factors. These 
solutions are known analytically and they are the 3 -I- I dimensional asymptotically AdS near horizon geometries [13] 


9nh 


dip'^ 

fW 





A = 


dt 
e — 
r 


( 8 ) 


where f{ip) ='0^ — 1—^ — ;| 2 , and where a, e and c are constants. We have taken A = —3. These solutions run 
from ipo <il} < oo, where is the largest zero of fipj), and the conformal boundary is at ip ^ oo, with metric 


9 = 


—dt^ + dr'^ 
^12 


a^dp)^. 


(9) 


The constant a is fixed in terms of e and c by requiring there to be no conical deficit in the bulk, and determines the 
conical deficit in the boundary metric. 

The special case tpn = a = 1 with e = 0 corresponds to pure AdS/i in the bulk, and the conformal boundary is 
asymptotically Minkowski with no conical deficit and no charge density. For the solutions we find here, we are either 
going to deform the boundary while keeping A = 0, or keep the boundary flat with no conical deficit while introducing 
an electric potential. We therefore elaborate on the extremal solutions in these two cases. 

For the case e = 0, we can write c = tpoitpo — 1). t/'o is the largest zero of f{'tp) so long as l/-\/3 < i/jq. The parameter 
a which keeps the bulk conical-singularity free can then be found to be a = ^'^^2 , so we have a one parameter family 
of solutions corresponding to the full range of possible values of . 







4 


If we instead impose a = 1, so the boundary is Minkowski, and allow a non-zero electric potential, we find another 
one parameter family of near horizon geometries. These solutions are discussed in m. where they are interpreted as 
point charge defects in the CFT. We find e in terms of the largest zero ipo to be 

e2 = V'^(l + 3^o)(l-V'o)- (10) 

In this case ipo is the greatest positive zero of / so long as 0.226 S 'tpo < 1. The parameter i/jq in this range labels 
a family of bulk solutions, and the dual CFT states are sourced by an electric potential e. As pointed out in 
the map, e —>■ ipQ, from source to state is neither left-complete nor one-to-one. The range of e for which there is a 
corresponding state is bounded by ^ (69 -I- ll-\/33) = e^, and, for ^ 0.066 = e^, there are two branches of 
states labelled by two different values of ■00 • The two branches are, 

• Branch 1: Runs from 0 < The = 0 case is pure AdS. 

• Branch 2: Runs from < e^. It meets Branch 1 at = e^. 

When we introduce a more general electric potential, we will again find two states for certain choices of sources. 


III. NUMERICAL SETUP 
A. Anzatz and Coordinate Fixing 

The task is to find 4 dimensional solutions of the Einstein-Maxwell equations, with conformal boundary given by 
([^, which are static and preserve the symmetry in ([^. These solutions generalise the near horizon geometries of (|^. 
We will take units where A = —3. 

The symmetries reduce this to an effectively two dimensional problem, with dependence on only an angular coor¬ 
dinate which we call 0 and a bulk coordinate which we call A, a compactified version of the tp coordinate from (|^. 
The dependence of the metric on the other two coordinates, r and t, are then fixed by the symmetries. We compactify 
the bulk coordinate, so that X lies in the range 0 < A < 1. A and (p taken together describe a unit disk, which we 
call M, with the boundary of the disk being the conformal boundary. 

For the bulk anzatz, we write down the most general metric and gauge field compatible with the symmetries 

1 / df‘^ dr^ dr ■\ 

9 = 7i-( 5'i(^,0)^ + S'2(A,0)— -h 2—uj{X,(j)) + g2ij{X,(j))dy^dy^ ) 

^ r r / 

A = 53(A,0)-, 
r 

where y = (A, 0). From the perspective of A4, Si, S 2 , and S 3 are scalars, w is a one-form, and g 2 is a symmetric rank 
2 tensor. Note that as we have explicitly factored out the divergent conformal factor, we can choose coordinates such 
that these are all smooth. 

We need to impose boundary conditions, and these need to fix the conformal boundary at A = I to be of the 
conformal class of ([^, and fix the electric potential to this value as well. To achieve this, the boundary conditions 
will satisfy 

.Sl(l,0)= 52(1,0) =47 
w(l,0) = 47(x(0)d0) 

92(1, 0) = 4(70^d0^ -f dA^) 

53(1,0) = R( 0 ), 

where we will choose the constant 7 as convenient. 

After fixing this anzatz we still have the residual coordinate freedom 

r ^ rA(A,r) 

A ^ A(A,y) 

Y Y{X,Y). 

We fix this using the Harmonic Einstein method. This method, and numerical approaches to applying it are discussed 
in [ini]. The essence is to gauge fix by choosing a gauge condition of the form 

e(a:)“ = g^-^ (F“^, - f= 0 , 


( 12 ) 


(13) 


(14) 
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where F is some chosen reference connection. This gauge condition is not imposed directly, but instead we add a 
term to the Einstein-Maxwell equations 




(15) 


= 0 . 


This additional term makes the Einstein equations elliptic for static spacetimes. In the absence of the gauge field, 
there is also a maximum principle that says so long as 4> = vanishes on the boundary it vanishes everywhere. 

This is derived in m by considering the contracted Bianchi identity. If the metric satisfies ( [l^ with F^i, = 0, this 
can be manipulated into the condition 


= -2A,^ + 2(V^r)(V^e.) 

In the static case, both terms on the right hand side are positive everywhere. This leads to the condition 


(16) 


(17) 


which, it can be shown, implies that (j) takes its maximum value on the boundary. Since, in addition, (j) > 0, this 


implies that if it vanishes on the boundary, it vanishes everywhere. From this it follows that if we solve (151 without 
a gauge field, with boundary conditions that ensure </> = 0 on the boundary, we should end up with solutions to 
Einstein’s equations. 

If F^i, ^ 0, we do not have such a maximum principle, and we will have to check that our gauge condition is indeed 
satisfied by the solutions we find to be sure that we are actually solving Einstein’s equations. Naively, you’d have to 
check for the vanishing of or equivalently (j>, everywhere. The issue with this is however small it is in the numerical 
solutions, you might worry that it doesn’t actually vanish completely. Fortunately, in some cases we can do something 


a little better. The addition of a stress tensor to Einstein’s equations, such as that due to F^v in (15), modifies the 


relation (16) to 


= -2A^ + 2(V'^r)(V^6) - 2T^.e'^e 


(18) 


where Note that in the static case, we can choose T such that ^ is spacelike, so that if the 

spacial part of —Agij — Tij is positive definite, the right hand side of this equation is again positive definite. We can’t 
in general assume that this will be true, but we can test whether this matrix is positive definite for the solution that 
we get, and, if it is, we know for sure that we must have ^ = 0. In principle, this can be a better test than directly 
testing (j) = 0, because, so long as the eigenvalues of the matrix don’t get too small, we can demonstrate that they 
are finite and positive without having to get into detailed arguments about convergence. 

In Figure we examine the positive definiteness for this matrix for the near horizon geometries. We see that 
throughout Branch 1, as defined in Section |II A[ this matrix is positive definite, whereas it becomes negative in 
Branch 2. For the more general numerical solutions, this picture isn’t changed too much, and there is a region of 
solutions where this matrix is positive definite, and a region where it isn’t. For the cases where this is positive definite, 
this demonstrates that we have actually solved Einstein’s equations, so long as (f) vanishes on the boundary, and, for 
the other cases, we have to check (j) directly. 

We will derive our reference connection T from a reference metric g which we will take to have the same general 
form as the metric So long as we fix the conformal boundaries of g and g to be the same, it follows, by examining 
the expansion of the condition (14) near X = 1, that we will have (^ = 0 on the boundary. We will therefore pick a 


reference metric which also satisfies (12). The precise choice will depend on the values of a, xW ^-nd V{(j)). 

Note that there is in general no guarantee that, for a given choice of g, the coordinate choice implied by the 
condition (14) is smooth. Indeed, as we show in Appendix in the case with a gauge field we can show that at the 
boundary our choice of coordinates is only smooth up to fourth derivative. That said, evidence from the analysis of 
the convergence of our solutions, presented in Appendix]^ suggests that our solutions in these coordinates are in fact 
C^. 


B. Discretization 


We have reduced the problem to a set of two dimensional pde’s, depending on X and </>, which should be thought 
of as polar coordinates on a unit disk with boundary conditions given by (12). Now we need to construct a lattice 
on which to discretize the problem. For simplicity, we would like to take a rectangular lattice, and we would like 
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FIG. 1: Minimum value of the minimum eigenvalue of the spatial matrix —A5* + Tj for the near-horizon geometries as a 
function of charge. For the entirety of Branch 1, as defined in this eigenvalue is positive, meaning that the matrix is 

positive definite. 


the conformal boundary, X = 1, to be one edge of this lattice. This leads us to using the polar X,(j) coordinates 
themselves as the lattice coordinates. 

These coordinates, of course, break down at the origin X = 0. In order to avoid this issue, we limit our attention 
to solutions that are invariant under parity 


(/)—>■(/) -h TT 


(19) 


To take advantage of this, we cover the disk twice by extending the X coordinate — 1 < X < 1. The extension is 
defined such that the map from polar to Cartesian coordinates is generalised to hold for negative X 


X = X cos 4> 
y = X sin (j). 


( 20 ) 


These coordinates map out the disk twice, and consistency then requires that when we extend our metric and gauge 
field to this double covering, they are invariant under 


(f) ^ (j) + IT, 


( 21 ) 


as this sends x —>■ x and y ^ y. Imposing this consistency relation would be messy, but, given our assumed parity 
symmetry (191, the solutions are in fact symmetric under each of X —> —X and (f> ^ (p + tt separately. 

Imposing each of these conditions separately is straight forward. We construct a lattice where X runs from —1 
to 1, with an even number of points so that we don’t need to place a point at X = 0, and p runs from 0 to tt. We 
construct a 6*’^ order finite difference matrix in the X direction, and then identify the points under X —^ —X so that 
we have resulting differentiation matrices for the points X > 0 which builds in the evenness without having a point 
at the origin which would require special treatment. In the p direction we use Fourier differencing to build in the tt 
periodicity. 
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In order for this to make sense, we need to choose an anzatz for the metric and gauge field in which the functions 
that define the components of both inherit the X —>■ —X and </)—)■ (^ + tt symmetries. Such an anzatz is given by 

S^{X,cj,) = {l + X^){l + A{X,cj,)) 

52(X,</>) = (1 + X2)(1 + B(X,</.)) 

uj{X, (j)) = F{X, (j))d{X^) + X^H{X, (j>)dcj) (22) 

g{X, (j)) = 4(1 + L{X, (j))){dx^ + dy^) + 2M{X, (j))dxdy + S{X, (j)){dx^ — dy^) 

S3iX,c^) = ViX,^), 

with X = X cos (j) and y = X sin (j>. In addition to the symmetry, this anzatz has two nice properties. Firstly, if all 
the functions vanish we are left with pure AdS. Secondly, a simple sufficient condition for these objects to be smooth 
on A4 is that the functions themselves should be smooth functions on A4. Note that, although we have avoided the 
numerical issues of having a point at the origin itself, X and 4>, being polar coordinates, are still not good coordinates 
at the origin. In order to test the smoothness of the functions at this point, we will need to transform back to the 
Cartesian coordinates x and y. This is elaborated on in Appendix [B| 

Having discretized the problem, we will find solutions using Newton itteration. An initial guess at a solution is 
made, and then we repeatedly improve our guess by solving linear systems. The linear systems are solved using 
UMFPACK[T 6 ]. 


C. Choice of Boundaries to consider 

Even after imposing (19), the space of possible boundary metrics and gauge fields of the form <§ is inhnite 
dimensional, so we need to focus on some particular subspace. We will separately consider the two cases where we 
deform the metric away from locally flat cones in the absence of a gauge field, and where we leave the boundary 
metric flat with no conical deficit, but introduce a non-isotropic electric potential. 

For the metric deformations, we choose the two dimensional subspace where 


X(((>) = Asin(2<))), 


(23) 


and we shall look at varying values of a and A. The one-dimensional sub-space, A = 0, is the set of boundary cones, 
for which, as described above, the bulks are the extremal solutions known analytically. Our reference metric needs to 
have this same deformed conformal boundary, so we take 

V , ^ ^ ^4.2 , , v2^m2\ , o \ v2 r, j.^rdcj) 


(i-x^y 


/ 1 (A)- 


(-C//2 + + 4 -H X^d(j)^) + 8XX^k{X) sin 2(j)- 


(24) 


with h{X) = 1-1- ^ X^ and k{X) = 1— (l — A^) . In order to impose the boundary condition that 5 ^ 1 / and g^i, 
should agree on the boundary, we impose the boundary condition with V{(jj) = 0 and 'y = a ^. 

For the gauge field deformations, we will fix a = 1 and x((/') = 0, and we will again take a two dimensional subspace, 
this time 


V{4>) = a + bcos{2(j)). 


(25) 


Again, the case 6 = 0 are the boundaries for which we know the two branches of analytic extremal horizon solutions. 
In this case the boundary is unperturbed, so we could in principle use pure AdS as our reference metric. However, 
we choose instead to use the two branches of extremal horizon solutions as reference metrics, as this makes it easier 
to find solutions close to each of these branches. We therefore write the reference metric as 


dfiv — 


{ 1 -x^y 




2(1 


A^ 


{-dy + dr^) + 4 + P(X)y^X^d(p^ 


(26) 


With PyX) = -- - - /2,y2_i_i,2 —^- -■ fhe parameter ipo is the parameter that labels these near 


horizon solutions from Section HA and it depends on the value of a through (10). In order to impose the boundary 
condition that g^^, and g^i, should agree on the boundary, we once again impose the boundary condition ( 12 ), this 


time with x(0) = 0 , a = 1 and 7 = ' 0 g. 

Note that, in addition to being invariant under ^ —>■(/)-|-tt as we required, the boundaries we’ve selected through our 
choice of s-nd Viy) have a (f> tt — (j) symmetry. Since we don’t break this symmetry in our choice of reference 
metric, and so in our choice of coordinates, it will be manifest in our solutions as well. 














AdS 2 Conformal Boundary 


r 


Axis X — 0 V 


Conformal Boundary X = 1 


Extremal Horizon 


FIG. 2: A Constant Time slice of the Bulk 


IV. RESULTS 

The results presented here have 25 points in the (j) direction and 145 or 65 in the radial direction for the boundary 
metric and electric potential deformations respectively"^. We already get good solutions at this resolution, but in 
Appendix]^ we demonstrate convergence by taking a subset of these solutions to higher resolution. For instance, in 
all regions of the solution space we explore, we are able in this way to push the error on the gauge condition, = 0, 
and on the Einstein’s Equations down to 0(10“®). This is discussed in detail in the Appendix. These solutions 
are found for a range of the parameters A and a in the case of boundary metric deformations, and a and b in the 
case of boundary electric potential deformations. They were found using the facilities of the Imperial College High 
Performance Computing Service [T7]. 


A. Electric Potential Deformations 

The boundary charge can be written as|llj 


Q = lim — 
z->-o 47r 




’S(z) 


(27) 


The memory requirements are higher with the addition of a gauge field. 
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where z is a function that goes to 0 on the boundary, and S(z) is a surface of constant z and t that becomes a spacelike 
hypersurface in the boundary as z —>■ 0. 

We’ve chosen coordinates where the conformal boundary at X = 1 is AdS 2 x but we need to be careful when 
relating this back to coordinates where the conformal boundary is Minkowski. When we consider things more carefully, 
we find that the full conformal boundary is divided into two components, one at X —>■ 1 which is the main component 
of the conformal boundary, and another at r —>■ oo, the conformal boundary of the AdS 2 factor in these coordinates, 
which is the piece of the boundary corresponding to the origin in the Minkowski frame. So, in the transformation 
from the Minkowski frame to the AdS 2 x frame, the origin has been blown up to the conformal boundary of the 

AdS2. 

A constant time, constant (p slice of our spacetimes is shown in Figure In this diagram, these two components 
of the conformal boundary are at the right hand side and the top. The full timeslice is the cylinder of revolution of 
this 2-d slice about the axis = 0. To calculate the charge as in ^7\ , we need to construct a sequence of surfaces 
S{z) that become a timeslice of the conformal boundary as z —>■ 0. Consider a cylinder, V, within the spaceslice as in 
Figure]^ where the top and bottom boundaries are at two values of r, ri and rg, labelled and respectively, 
and the other surface of the cylinder. Ex, is at a fixed value of X. In the limit that we take this cylinder large. 
Ex U Er^ approaches both components of the conformal boundary. It is therefore a suitable choice to extract the 
boundary charge. We can write 


Q= lim f *F+ f I . (28) 

X-)-l,ri-i.oo 47r Jsx ) 

These two terms represent the contribution of the two components of the conformal boundary to the total charge. 
From the point of view of the Minkowski frame, the first term is a delta function charge density at the origin and the 
second term represents a charge distribution away from the origin. 

In fact, there is no overall contribution from this second term. We see this as follows. The Maxwell Equation can 
be written 


d*F = 0. 

We can integrate this over the volume V using Stoke’s theorem. 


(29) 


d*F = 


*F + 


*F + 


*F = 0, 


where all the surfaces are oriented outward from V. 


Our scaling symmetry in fact implies that 




*F, 


so that. 


(30) 


(31) 


[ *F = 0. 


(32) 


This means that the only overall charge contribution can come from the origin, although there will be a non-zero 
charge density elsewhere on the boundary. In Figure]^ we plot some examples of charges at the origin as a function 
of a for various values of b for the boundary source potential V{4>) = a + bcos2(j), together with the analytic curve 
for the near horizon solutions (6 = 0). These curves all go through the origin, which makes sense since, when a = 0, 
there’s nothing to pick out a preferred sign of charge. For each value of a and 6 we found these solutions by using 
a corresponding near-horizon solution with the matching value of a as the reference metric and as our initial guess 
for the metric, and then proceeding with a simple Newton itteration. In the regions of values of a where there are 
two branches of near-horizon solutions, we tried using each of these branches. In this way, we were able to find two 
solutions for a range of a and 6 pairs. For instance, we can see in the b = 0.8 curve in Figure]^ a range of a values for 
which, like in the near-horizon case, there are two branches of solutions. For the larger values of 6 we weren’t able to 
find two branches of solution in this way, but we have no reason to rule them out. Additionally, we weren’t able to 
find solutions for values of a greater than the maximum value on the near-horizon solutions. 

In addition to the delta function charge, there is a non-zero charge density away from the origin, which in our case 
corresponds to a non-zero charge density on the X = 1 conformal boundary. In the AdS 2 x frame, because of the 
scaling symmetry, both the electric potential and charge density fall off as ^ , and have non-trivial dependence only 
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- b = 0.8 

. b = 1.6 

- b=2 

. b=3 

- near horizon analytic solutions 


0.0 0.1 0.2 0.3 0.4 0.5 0.6 

a 

FIG. 3: Overall charge on solutions as a function of the parameters a for a section of choices of b in the electric potential source 
term V{(j>) = a + b cos 2(j> for which we find solutions. Note that at a = 0 the total charge is always zero, which we would expect 
from symmetry. 



on (j). In the Minkowski frame this corresponds to a l/r fall off for the potential, and a 1/r^ fall off of the charge 
density. From (32) we expect the integral of this charge density to vanish, and in all of our numerical solutions this 
integral does not rise above 10“®. In this large scale region we also look at the induced energy density. This energy 
density is found via a Fefferman-Graham expansion [181 [19], and, in the Minkowski frame, it falls off as 1/r^. 

We show plots from two sets of examples to illustrate how these field theory observables depend on the source 
potential. First, the a = 0 case where V{(j)) = & cos 20 for various values of b . This is shown in Figure]^ In this case 
there is only one branch of solutions, and there is a symmetry between the regions of positive and negative charge 
density because our boundary conditions don’t pick out a sign of charge. The energy density is positive definite and, 
as you might expect intuitively, it’s largest in regions of large density. If we consider instead F(0) = 0.5+ 6 cos 20, in 
Figure]^ we find two branches of solutions. This time, because of the offset of F(0) from zero, the symmetry between 
regions of positive and negative charge is broken. This can be seen most clearly in the energy density distribution 
which is larger in the region of positive charge density than in the region of negative density. Once again, the energy 
density is positive definite, and in fact seems to stay above the pair of constant values it takes in the corresponding 
6 = 0 extremal solutions. 


B. Metric Deformations 


The other set of solutions we found were ones where we did not introduce a gauge field, and only deformed the 
boundary metric. In this case, the observable we focus on is the X = 1 boundary stress-tensor, which we extract 
using a Fefferman-Graham expansion. We will focus on the energy density, which, as remarked above, goes as I/r^. 
Whereas in the case of the electric potential deformations the energy densities were positive, once we deform the 
boundary metric we get regions of negative energy density. This was already true for the near horizon solutions, 
described in Section [ll A All we do in that case is introduce a conical deficit on the boundary (x(0) = 0). The 
resulting Ttt{4>) is a constant which is positive for a > 0 and negative for a < 0. 

Once we break the iS'0(2,1) symmetry by introducing x(</’) = A sin 20, this becomes non-isotropic and there can 
be regions of positive and regions of negative energy density. In Figure we give some examples of energy density 
profiles for a given set of profiles. The energy densities are shown at three different values of a, one corresponding 
to no-conical deficit (a = I), one to positive conical deficit, and one to negative conical deficit. Qualitatively, the 
behaviour doesn’t change too much between the three cases. The energy density we get oscillates about the ^(0) = 0, 
near-horizon, value with an amplitude that grows as ^ grows. The offset from the near-horizon value is largest when 
the gradient of ^'(0) is largest, and vanishes close to the point where C^(0) = 0. 
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(a)Form of the applied source V{ip) 


(b)Resulting charge density 



(c)Resulting energy density 


FIG. 4: Response to an electric potential source term for b — {0, 0.8,1.6, 2.4, 3.2}. The induced charge density 

integrates to zero, but the energy density is positive definite. Note that, quite intuitively, the energy density is peaked when the 
the magnitude of the charge density is largest. Also, there is a symmetry in the energy density between the regions of positive 
and negative charge, which makes sense because there is no overall charge at the origin. Viewed in a flat frame the potential 
goes as 1/r, the charge density as 1/r^, and the energy density falls off as 1/r^, and we have evaluated all of them at r = 1. 
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(a)Form of the applied source V((/>) 



(d)Resulting energy density Ttt{(j>) for branch 1 (e)Resulting energy density Ttt{4>) for branch 2 


FIG. 5: Response to an electric potential source term V{(j)) = for h = {0, 0.2, 0.4, 0.6,0.8}. In this case we found 

two branches of solutions. Note that even though the potential V{^) is no longer centred about zero, unlike in Figure l^the 
charge density still is, and it still integrates to zero. The main difference now is that we no longer have a symmetry in the 
energy density between regions of positive charge and regions of negative charge. This is because these solutions have an overall 
(positive) charge at the origin. Again, they have all been evaluated at r = 1. 
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E E in E. 

8 4 8 2 


(a)Form of the off-diagonal SO(2,l) breaking drdcf) term in the boundary metric 

x{4>) 



(b)Resulting energy density Ttt{(i>) for a = 1 




(c)Resulting energy density Ttt{<f>) for a = ^ 


<P 

(d)Resulting energy density Ttt{(p) for a 


fs 

V 6 


FIG. 6: Large scale energy density for a boundary metric with x{4^) — ^ ~ 0.01,0.05, 0.1, 0.15}, with 3 different 

conical deficits a = energy densities go as 1/r^, and we’ve evaluated them at r = 1. Now that we’ve 

deformed the boundary metric, the energy density is no longer positive definite. This is not a feature that is peculiar to our 
solutions as it was already the case in the near horizon solutions that some had positive energy density, and others negative. 
Roughly speaking the energy density follows the gradient of x(0)? i't always intersects the near-horizon energy density for 
the corresponding value of a close to when — 0- 
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V. THE SINGULARITY 


To understand the nature of the singularity that replaces the extremal horizon, we examine uncharged geodesics 
that approach this surface. Moving, for the moment, to general dimensions d, we choose coordinates where the metric 
takes the form 


_—I- . dr 

ds^ = -s-f 2'>p{x)‘^Ai{x)dx^ -h {hij{x) + 'ip{x)‘^Ai{x)Aj{x))dx''dx^, 


(33) 


so that the null surface of interest is now at r = 0. The a;* coordinates describe a d — 2 dimensional subspace, on 
which hij(x) is a non-degenerate metric, and Ai(x) and &re a one-form and scalar respectively. The near horizon 
solutions have Ai(x) = 0 when written in this form, and the generalized solutions we found above have no-zero Ai(x) 
because they break the SO(2,l) symmetry. Using the remaining global scaling and time translation symmetries, we 
can integrate some of the geodesic equations using Noether’s theorem. The equations we are left with are. 


E = 

K = Et + \ + Aix‘ 


(34) 


= - ^ (fc2 -K hkix'^x^) + EAH - E)x^ {K - Et) 


where a dot represents differentiation with respect to an affine parameter A. Indices are raised and lowered using hij, 
and we have defined Fij = Aij — Aj^i. K and E are conserved quantities due to the static and scaling symmetry, and 
the third equation is the geodesic equation in the a;* direction. The parameter is 0 for null geodesics, and can be 
set to 1 for timelike and —I for spacelike geodesics. 

By dotting the third equation in with we can show that® 


dX 


+ XjX^)) = 


2E^x^Ai 


Substituting this equation into the conservation equations, we find the differential equation 

1 d 


which can be solved to give 


K \ dr 

E ^ dt 


\t) = 


2 U2 dt 


-I- XjX^)) , 


E^ 


{Et -K)^ + C + V'2(fc2 -h XjX^i) ■ 


(35) 


(36) 


(37) 


The parameter C is a constant of integration. From this equation it can be seen that any time-like or null geodesic 
will reach r = 0 as t —)■ oo. Using this result, we can rewrite the third equation as 


= (^A^ - + hkix'^x‘) + ^ {{Et - Kf + C) + F^x^ {Et - K). 


(38) 


Now, assume there exists some point on the n — 2 dimensional x manifold, xg, where Ai = dipi = 0. Such a 
point exists, for instance, at the point A = 0 in the solutions we constructed, because of the reflection symmetry we 
imposed. We can then consider freely falling particles at this point with i® = 0. In this case, the particles stay at 
X = Xg. Timelike geodesics satisfy 


,,, cot A 
t{X) =-——I- to 

r(A) = —BsinA, 


(39) 


This becomes an extra conservation law when Ai = 0, which is because the scaling/static symmetry gets enhanced to S'0(2,1) in this 


case. 
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Iq and B being constant, while null geodesics are given by 


t{X) ~ ~ X to 
r{X) = -A, 


( 40 ) 


where in both cases we have chosen A = 0 to be the time when the geodesic hits the surface. 

Define 14“ to be the tangent vector . The geodesic deviation tensor in the n — 2 dimensional x* space is then 
given by Rij = RaibjV°‘V^. This means that if the vector X* describes the infinitesimal displacement between a pair 
of geodesics both initially parallel to 


V.V{V.VX^) = -h^mjkX^. 


(41) 


For both timelike and null cases we find that this goes as 


R 




(42) 


We therefore see that, at least for this geodesic, at ^^(x) = d^{x) = 0, the geodesic deviation diverges as you approach 
the surface. 

More generally, the geodesics seem to be very poorly behaved. As r —)■ 0, we have t —>■ oo. By examining (381, 
it seems to be difficult to have this, and not have divergent x, unless A* = 0 at the point the geodesic hits r = 0, 
although we note that we don’t have a proof of this. 

We can, however, provide evidence for this numerically. We can take our numerically found bulk geometries (with 
no gauge-field), and shoot timelike geodesics towards the singularity. In Figurej^we plot x(r), and x'(r) for a selection 
of geodesics on pure-AdS and on one of our 5*0(2,1) symmetry breaking solutions. The coordinate x is one of a pair 
of Cartesian coordinates on our unit disk, and r is proper time. These geodesics were evolved using an adaptive step 
size Runge-Kutta method, and we integrated until the step-size became too small to continue. From this we can see 
that, in the non-AdS case, as the geodesics get close to the singular surface there is a large transverse acceleration, 
and it looks like geodesics will have divergent x by the time they hit the surface. It seems reasonable to conclude, 
based on examination of (38) above, that the large transverse acceleration observed in the twisted case just before 
the end will continue and 


ead to a divergent x. 




16 



(suJAdSi 



(b)SO(2,l) breaking solution o = 1, A = 0.15 



- 0.020 - 0.015 - 0.010 - 0.005 0.000 


{c)AdS 4 , 



(d)SO(2,l) breaking solution 


FIG. 7: x{t) and x'{t) for various numerical timelike geodesics that head towards the r — 0 surface. The proper time coordinate 
has been chosen so that it is 0 when they hit the surface. As they approach the surface, nothing happens in the case of the AdSi 
and they just cross over a smooth extremal horizon. In the other case, all the geodesics start to swerve off as they approach 
the horizon, except a single geodesic that we’ve placed at the point x = y = 0 with no transverse velocity (*' = t/' = 0). This 
is the point where Ai = dtpi = 0 as discussed in the text. 


VI. SUMMARY 

We have found 3+1 dimensional static, scale-invariant solutions to Einstein-Maxwell equations outside the class of 
the known analytic solutions which correspond to the possible near horizon geometries of extremal horizons. We would 
argue that these are the appropriate geometries to describe the infra-red when we are considering states in holographic 
field theories on 2+1 dimensional boundary geometries, which, on a large scale, do not have the full SO(2,l) isometry 
needed to produce an extremal horizon in the bulk. We would expect that if you took these non-scale-invariant 
geometries, and took a large scale limit, you would end up with the scale-invariant geometries we found here. 

For instance, we can consider introducing a source electric potential into a holographic CFT which has some arbitrary 
form, but falls off as A ^ As remarked in m, we can use extremal horizons to describe such electric-potentials 
which asymptotically fall off as 


A = 


edt 

r 


However, if the fall off breaks this axial-symmetry. 


_ V{0)dt 

r 


(43) 


(44) 


the SO{2, 1) symmetry is broken, and we expect the singular solutions we’ve found to be the relevant infra-red 
geometries. In particular, we would expect there to be field theory states where the expectation value of the charge 
density falls off as p{(j))/r‘^, where the p{4>) has the form determined by the scale-invariant geometries we’ve found here. 
Similarly we’d expect them to have an energy density that has a 1/r^ fall-off again determined by these scale-invariant 
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geometries. Note that, for the near-horizon solutions, there is no such boundary charge density away from the origin, 

p{4>) = 0. 

These charge densities and energy densities are the ones we gave examples of in Figures and (and also in Figure 
l^for the case of metric deformations). We also looked at a delta function charge at the origin of our scale-invariant 
states in Figure]^ We would expect this in usual circumstances to correspond to the overall charge of the non-scale- 
invariant states of which these are the large scale limit. Indeed, it is shown in that if you start from no chemical 
potential in a CFT, and smoothly turn a chemical potential on with an arbitrary dependence on r, then it is only 
a marginal chemical potential which leads to a finite net contribution to the total charge. This marginal chemical 
potential is the scale-invariant A ~ Vdt/r which we have in our singular geometries. 

One limitation we ran into was the lack of smoothness in our coordinate choice at the boundary. This is explained 
in Appendix]^ and it stems from our choice of reference metric. An interesting question for future work using the 
technique of the harmonic Einstein equations would be how to pick a reference metric that leads to coordinates that 
are as smooth as possible. 


Acknowledgements 

We would like to thank Toby Wiseman for many helpful discussions, advice, and feedback. The numerical work was 
carried out using the facilities of the Imperial College High Performance Computing Service. All plots were produced 
using Mathematica. AH is supported by an STFC studentship. 


Appendix A: Gauge Fixing Static Scale Invariant Geometries 

Here we explain why the most generic 2-1-1 dimensional static, scale-invariant boundaries, where the scaling 
symmetry takes the form 


t —^ ct 
r —>■ er, 


(Al) 


can be written in the form @. 

First we consider the metric. If we write down all the terms in the metric that satisfy the symmetries, we get 


fl'p 

g = —-I- -\- 2H{(j))dip —. 


(A2) 


We are left with the following coordinate freedom 

i Weyl Scaling of the metric since we only care about it’s conformal class: g' = exTp{f{(j)))g 

ii The redefinition of the r coordinate: r = X{4>)r' 
hi The redefinition of (p: (p = (p{(p'). 


Note, that there is an addi tiona l freedom corresponding to the rescaling t = /it'. However, because of the symmetry 
of (Al), this is equivalent to|(ii) with X{p) = i. 

and (ii) to set A = H = 1. This uses up these two coordinate freedoms completely. 


We can use the freedom of 
to set 


We can then use (hi) 


ad(p' = S{p)dp. 


(A3) 


We choose to fix cp and p' to have canonical period 27r. We are therefore left over with the constant cn = -^ J S{p)dp 
because we are restricted to redehnitions of p that preserve this period. 

Dropping the prime we hnd that we can write the most general boundary as 


dv dr 


dr 


g = -y H-^ -I- a^dp^ + 2ax{p)dp —. 


(A4) 


Therefore, a 2-1-1 dimensional twisted cone boundary is specihed by one periodic function x(0) and one constant a 
which we take to be positive. 
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If we demand that these coordinates are non-degenerate, we require 


lx(</>)l < 1, 


(AS) 


so that det( 5 ) < 0. 

To check that these boundaries aren’t secretly conformally flat, we look at the Cotton-York tensor. The non¬ 
vanishing components of this symmetric, traceless tensor are 


^ 2x(</>)x/(0)^ _ X"i'p) 

a(x(0)^-l)^ 2a (x((5!))2 - 1)^ 

ar{x{4>y - 1)^ 


(A6) 


The vanishing of this tensor is the necessary and sufficient condition for conformal flatness in three dimensions. We 
see that for the case of constant x{4‘) = M this condition is satisfied, and there must be some extra gauge freedom. In 
this special case, we can make the redefinition (j) ^ (j) + /rlogr to get rid of x(<^)^- 

The gauge field is simpler. We have a t ^ —t symmetry which means that it is consistent to consider a purely 
electric source potential, with no spatial vector potential. When combined with the scaling symmetry, this allows us 
to immediately write down 


A = V{ct>)-. (A7) 

r 

Note that we cannot use the redefinition of r, |(ii)[ to get rid of this, as we have already used up this freedom in fixing 
the metric. 


Appendix B: Convergence 


In order to make the claim that the solutions we find are good approximations for actual continuum solutions, 
we need to have some checks of convergence. We perform these checks using 6**' order finite difference in the radial 
direction, and spectral interpolation in the angular direction, for a selection of the a, b pairs we collected for the 
solutions with a gauge field. We also performed corresponding checks for the pure gravitational solutions and we see 
similar behaviour. 

We start by checking convergence of our solutions directly. Plotted in Figure]^ is the integrated absolute value of 
the shift of one of the functions in our solutions between successive X resolutions at a fixed (p resolution of 20 points. 
Comparing these to the corresponding plots in Figurej^ where we fix the X resolution to 65 and vary the (p resolution, 
we see that’s its the dependence on X resolution that dominates as the shifts in this case are orders of magnitude 
bigger. 

We therefore focus on dependence X resolution. Since we are using 6*^ order finite difference, we might expect 
these shifts to fall off as (so that the offset from the continuum values falls as n“®), where n is the number of 
points. However, what we instead see is an n~^ fall off in the shifts of the function and its first derivative, and then 
more like n~‘^ for it’s second derivative. This is consistent with the limitation of smoothness at the boundary coming 
from our choice of coordinates, as elaborated in Appendix [E} Indeed, in Figure [T0| we plot shifts in just the boundary 
values of the field, and we see slower convergence. The fifth derivative is not converging at all, and it is at this order 
that the boundary expansion in Appendix [E| tells us the lack of smoothness gets introduced. 

As a further check, we can examine our gauge condition. The vector = g°‘^ (~ ) should be converging 


to zero. As in our case it’s a space-like vector, it is sufficient to look at p = It’s dependence on both X 

resolution and p resolution is shown in Figure El These seem to demonstrate convergence, although in the latter 
part of 11(a) in which we plot dependence on X resolution, we are limited by the fixed p resolution, and vice versa 
for |ll(b) Since our solutions satisfy the harmonic Einstein equation, this is an indirect test of how well they satisfy 


the Einstein equations. We see the same behaviour when we examine the trace of the Einstein equation. 


As is evident from Equation |A6| we cannot in general use a redefinition of this form to get rid of x- If we set rp (pX xi'P) fog r we pick 
up log terms in the metric unless xi'P) constant. 
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• a = 0., b = 2. 

■ a = 0.2, b = 1.5 

♦ a = 0.4, b = 1. 

A a = 0.6, b = 0.5 

V a = 0.65, b = 0.3 



— 1.16795 -5.25177 X 

— 0.994211 - 5.29023 X 
■— 0.815286-5.32267 X 

— 0.792797-5.41743 X 

— 0.972055-5.56309 X 


® a = 0., b = 2. 

■ a = 0.2, b = 1.5 
♦ a = 0.4, b = 1. 

A a = 0.6, b = 0.5 
T a = 0.65, b = 0.3 



1.339 -5.24203 X 
1.22197 -5.30239 X 
1.03653 -5.33225 X 
1.09134 -5.44232 X 
1.46479 -5.65005 X 


. a = 0., b = 2. 

■ a = 0.2, b = 1.5 
A a = 0.4, b = 1. 

A a = 0.6, b = 0.5 
T a = 0.65, b = 0.3 


Logio IntegratedDifferencein 



Resolution 


1.5719 -4.26589 X 
1.71444 -4.39728 X 
1.67734 -4.46087 X 
1.37568 -4.43838 X 
1.15202-4.42301 x 


FIG. 8: Integrated absolute value of the shift in the solutions between subsequent X resolutions with 20 points in the (j) 
direction, taken on the slice (p = ^. We take the case of the function A and we plot the shifts in (a) it’s value, (b) it’s first 
derivative, and (c) it’s second derivative. These are log-log plots, and our linear fit to the final few points give us an estimate of 
the order of convergence. In these fits, x = log^gResolution. These seem to indicate fourth order convergence for the function 
and it’s first derivative, and second order convergence for it’s second derivative. 
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Logio Integrated Difference in A 

; 5 10 

-5 

ff 

• a = 0., b = 2. ♦ • 

■ a = 0.2, b = 1.5 * ♦ 

T 

♦ a = 0.4, b = 1. ,1 

» a = 0.6, b= 0.5 ^ 

" a = 0.65, b = 0.3 

-25 

-30 

(a) 


Resolution 

15 20 25 


Log-jo Integrated Difference in A' 


• a = 0., b = 2. 

■ a = 0.2, b = 1.5 

* a = 0.4, b = 1. 

A a = 0.6, b = 0.5 
T a = 0.65, b = 0.3 


(b) 


Logio Integrated Difference in A” 

5 10 

-5 . 

♦ • 

. a = 0., b = 2. i 

-10 . «. 

■ a = 0.2, b = 1.5 

♦ a = 0.4, b = 1. * 

i a = 0.6, b = 0.5 

V a = 0.65, b = 0.3 

- 20 - 

- 25 - 

(c) 


' Resolution 

15 20 25 


I 


FIG. 9: Integrated absolute value of the shift in the solutions between subsequent 0 resolutions, with 65 points in the X 
direction, taken on the slice (f) = We take the case of the function A and we plot the shifts in (a) it’s value, (b) it’s first 
derivative, and (c) it’s second derivative. These plots indicate approximately exponential convergence with (j> resolution. When 
you continue above a 4> resolution of 25 the X resolution becomes limiting. 
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• a = 0., b = 2. 

■ a = 0.2, b = 1.5 

♦ a = 0.4. b = 1. 

i a = 0.6, b = 0.5 
a = 0.65, b = 0.3 


Logic A'[ 1] 



1.00518 - 4.07895 X 

— 1.30525 -4.28402 X 

— 1.32121 - 4.37266X 
- 1.14994-4.40247 X 

0.963164-4.40154X 


Logic A"[1] 



1.41634 -3.06403 X 
1.80379 -3.30672 X 
1.85182-3.40881 x 
1.6927 -3.44383 X 
1.50892 -3.4442 X 


• a = 0., b = 2. 

■ a = 0.2, b = 1.5 

♦ a = 0.4, b = 1. 

i a = 0.6, b = 0.5 
T a = 0.65, b = 0.3 


Logic/A(^)(1) 



—■ 1.95268 -0.188868 X 

— 2.44509 - 0.47768 x 

— 2.53009 -0.595595 X 

— 2.39088-0.639306 X 

— 2.22168-0.646226 X 


FIG. 10: Absolute value shift in the boundary value of the solutions between subsequent X resolutions, taken at (j> = 

We take the case of the function A and we plot the shifts in (a) it’s first derivative, (b) it’s second derivative, and (c) it’s 
fifth derivatives. Linear fits are made to the final few points with x = logj^gResolution. Convergence is slower than the bulk 
convergence in Figure and at fifth order we see the effect of the lack of smoothness explained in Appendix [E] 
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FIG. 11: The maximum value of <j) = S'S a function of resolution. In (a) we see the dependence on X resolution at a fixed 
0 resolution of 20 points. Initially we see what looks like power law behavionr in all these solutions, but a couple of them 
start to level off at around 65 points. This seems to be where the 0 resolution is becoming the limiting factor. In (b) we plot 
the dependence on (j) resolution at a fixed X resolution of 65 points. Above the (j> resolution of 20 which we plot up to the 
X resolution becomes the limiting factor. For the parts of (a) which look like a power law, the slope is about —3.9, which is 
consistent with the third order convergence seen in the second derivatives in Figure]^ 
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We have discussed the limitations of smoothness at the boundary, but we also need to check how smooth our fields 
are in the bulk^. The convergence plots in Figure give us some degree of confidence that the metric functions are 
at least C'^, but the coordinates that analysis is done with respect to are polar coordinates, so this might go wrong 
at the origin. To check against this, we move to coordinates that aren’t degenerate at this point. 

The simplest choice of such coordinates are the Cartesian coordinates {x,y) related to these polar coordinates in 
the usual way. We therefore transform the scalars iSi, S 2 and V, the vector oj, and the metric 52 as defined in (111, 

into these coordinates, and examine an expansion of their components near the origin x = y = 0. 

If a function, such as the components of these objects, is a smooth function of x and y at the origin, this translates 

into a particular dependence on X and cj) there. We can expand such a function in Fourier modes 


fiXcj}) =^/„(X, (ji) = (a„(X) cos ncj) + bn{X) sinnt/)), 


(Bl) 


and see what happens when we impose smoothness mode by mode. The symmetries X —>■ —X and 0 + tt of our 
solutions simplify things, and mean that the scalar and tensor components consist only of even modes of the form 


/e(X,</>)=^/2„(W</')=E {A 2 n cos2n^ + B 2 n sin 2n(/)) 


(B2) 


whereas the vector components are the odd modes 

fo (W, 4>) = f^n+i [X, 4>) = X (^ 271+1 {X'^) cos(2n + l)(j) + B2n+i (X^) sin(2n + 1)(^) . (B3) 


Note in particular that, thanks to the symmetry, the functions Aj and Bj depend only on X^. Since a Fourier mode 
like this is a smooth functions of the related Cartesian components x and y if it takes the form 

fn {X, (j)) = X'^ (^An (X^) cos net) + Bn (X^) sin nej^ , (B4) 

all we need to check for each Fourier mode is that they fall off as X”. We can turn that on it’s head and demand that 
if we expand our functions in powers of X, the (j) dependence at order X" should be a combination of terms of the 
form cos mej), and sinmi^ for m < n. If we show that this holds up to order X"*, we will have shown that the functions 
have well defined derivatives in x and y up to order m at the origin. In Figure we check for Fourier modes that 
violate this condition at order X°, X, and X^. We see very fast convergence in the first two plots, and even at order 
X4 we are still seeing fourth order convergence. This seems to indicate that the functions are well behaved at least 
up to the fourth derivatives at the origin. 

Another check of agreement with Einstein’s Equations is to check the asymptotic expansion of these solutions near 
the conformal boundary. To this end, we ch eck whether the extracted stress tensor is traceless and satisfies the correct 
conservation law®. This is shown in Figure [T^ We see that these quantities do not converge very rapidly, but, given 
the lack of smoothness on the boundary, we would not expect them to. 


^ As described in Section this manifold is singular, however we can still impose smoothness away from the singularity r = 0. This 
amounts to demanding that, as defined in the scalars and S 2 , the vector a;, and the metric g 2 are smooth over the 2-d disk A4 

® The presence of the gauge field source term in the CFT means that the conservation law for the stress tensor becomes 
where the factor | just comes from the fact that we have used a non-standard normalization of the stress tensor. 
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• a = 0., b = 2. 

■ a = 0.2, b = 1.5 
. a = 0.4, b = 1. 
i a = 0.6, b = 0.5 
a = 0.65, b = 0.3 



3.55104 -7.95254 X 
3.02816 -7.9929 X 
2.59547 -8.02795 X 
2.47398 -8.07289 X 
2.39167 -8.04304 X 


^ Logio 1 1 mode] ^ 
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• a = 0., b = 2. 

■ a = 0.2, b = 1.5 

♦ a = 0.4, b = 1. 

i a = 0.6, b = 0.5 
T a = 0.65, b = 0.3 



— 3.15515 -7.921 x 

- 2.68844 -7.94865 X 

■— 2.22616 -8.0107 X 

— 1.72826 -7.92311 X 

— 0.649843-7.54824 X 



4.82329 -3.92457 X 
4.20472 -3.93715 X 
3.39765 -3.94744 X 
2.44707 -3.97635 X 
1.82727 -4.07136 X 


FIG. 12: We’ve expanded our fields near the origin X = 0 in a Cartesian basis {x,y) and checked that derivatives in x and y 
are well defined. Thanks to our parity symmetry, this amounts simply to checking that at each order X”, if we expand the (j) 
dependence in Fourier modes, there are no modes above cos n(f) and sinm^. We plot here the sum of the squares of the modes 
cosmc^ and sin m(j> at order X" for m > n summed over all the components of our fields, for (a) n = Q, (b) n = 1 and (c) n = 4, 
at a series of X resolutions, at a fixed 4> resolution of 15 points. These show that the functions are well behaved at least up to 
fourth derivatives at the origin. 
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• a = 0., b = 2. 

■ a = 0.2, b = 1.5 
a = 0.4, b = 1. 

A a = 0.6, b = 0.5 
T a = 0.65, b = 0.3 
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(a)Tracelessness 



2.49037-1.91035 X 
2.2169-1.88986 X 
1.91577-1.87052 X 
1.53822-1.8546 X 
1.24973-1.8422 X 


FIG. 13: We plot how well the stress tensor satisfies the conservation and trace requirements. Neither converges very rapidly, 
but given the lack of smoothness on the boundary, we would not expect them to. 


Appendix C: Check of numerics against Extremal Solutions 


As discussed in Section II A[ when xW = 0 and the gauge field is rotationally symmetric, the static bulks that 
preserve the symmetries of Equationare smooth 4 dimensional near-horizon geometries. These solutions are known 
analytically, with their form given in Equation We focus only on the pure gravity case here, so no gauge field, and 
use these solutions as a sanity check. 

In this case, Weyl squared is given by 


2 _ 12 -'ipoY 


C" = 


'(/j® 


(Cl) 


In addition, the grr component of the metric is invariant under the coordinate freedoms 13 so we know that 

= (C2) 
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We can combine these together to get a scalar relation which should be satisfied by our solution 

C^gl = 12(^1 (C3) 

The largest error we find on this comes near the axis of symmetry, and we plot a comparison for a couple of choices 
of a in Figure [Til 


Appendix D: Comparison with Linearised Solutions 

In (rUj we presented a linearised calculation for the case of no gauge fields. We reproduce it here, with the addition 
of a gauge field, for completeness, and so that we can compare it to the numerical solutions we have presented. We 
look at linear perturbations of AdSi that preserve the static scaling symmetry. Because of the rotational symmetry 
of the background we’re expanding about, we can consider Fourier modes and write a gauge fixed perturbation in the 
form 
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9=^ {-dt^ + dp^) + “ l) 

+ ecosn(/i ^p„(V') + {'•P^ - 1) dtj)^ 

+ e cos ncj) f 2 /„— 

\ P 

+ esin710 ( 2/„('(/))—c?0 
V P 

+ e sin 710 (2/i„ {'ip)dij^d4>) 

A = €Vn{tp) cos 770 — 

P 

where n > 0. The gauge field decouples and just satisfies Maxwell’s equation on the ^^54 background. The solution 
is 




VnW = Bi 


Ip — 1 
Ip + 1 


njl 


Bo 


Ip — 1 
0 + 1 


—nj2 


(D2) 


where Bi and B 2 are constants. This has a divergence at the origin unless B 2 = 0. This gauge field makes no linear 
order contribution to the stress tensor, but it leads to a charge current on the boundary 


J^(0) = —enBi cos 770— 


(D3) 


For the gravitational perturbation, we expand Einstein’s equations to linear order, the rest of the problem reduces 
to a single second order ode for p„, with the other functions determined in terms of this, with solution 


Pnilp) = 


1 


0(02 _ 1 ) 


Cl 


Ip — 1 

0+1 


{nip -\) + C 2 


0+1 

0—1 


{nip + 1) 


Requiring that the perturbation is smooth forces C 2 = 0. The other functions are then 

ln{lp) = -Cinip{ip - 1)5“\0 + 


/n(0) = 
/l„(0) = Cl 

The boundary metric gets deformed to 


0—1 


n/2 


(nip + 0^ — 1) 


0 + 1 , 

(0 — l)t“i(0 + l)“t“i(n,'0 — 1) 

0 


ga = 


-dt^ + dp^ , ,2 „ „ . , dp ,, 

- + d(p^ + 2eCi smncp—dcp, 


P P 

and, through the Fefferman-Graham expansion [19). we can extract the boundary stress tensor 


0 0 


^0 0 0 


= C'iecos770 


" 0 + Cle sin 770 0 0 ^P] ■ 


0 -|y 


,0 0 

‘ 3p 


(D4) 


(D5) 


(D6) 


(D7) 


We compare these charge and energy densities to the ones we found numerically in Figures [Tb] and [Tfil As one would 
expect, initially they agree well, but there is a disagreement that grows once the perturbation grows to big. 
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<P 


Linearized Solution 
Numerical Solution 


FIG. 15: Charge density with an electric potential source term V (</>) = ^ for a range of values of b compared to that of 
the linearized solution 



0 


Linearized Solution 
Numerical Solution 


FIG. 16: Energy density for a boundary metric with x{4') 


^ compared with the linearized solution. 


Appendix E: Logarithm in the Boundary Expansion 

We derive here a limitation on the smoothness of our choice of coordinates at the conformal boundary for the case 
of our solutions in the presence of a gauge field. We start by assuming that the metric is smooth when written in 
Fefferman-Graham coordinates. In these coordinates, asymptotically AdS solutions to the 3 + 1 dimensional Einstein’s 
equations can be expanded near the conformal boundary as [19] 

dz^ + {gij{x) + rij{x)z^ + Tij{x)z^ + O ( 2 ^)) dx^'dx^ 

■ (El) 


The conformal boundary is at z = 0, gij is the metric on the conformal boundary, is determined by gij, but 
Tij can’t be fully determined by an asymptotic expansion, and is proportional to the induced stress tensor on the 
boundary. 
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Our solutions are not written in these coordinates, but instead in coordinates fixed by the gauge condition 

e = (E2) 

The reference connection, is derived from a reference metric gfj,u, which we fix in a advance. 

In our case, we have chosen to itself be a solution of Einstein’s equations, which means that it can also be 
written in Fefferman-Graham coordinates, in which it has the expansion 


dz'^ + i^gijix) + rij{x)z'^ + Tij{x)z^ + O {z^)) dx''dx^ 
= ^2 

In order to impose 0 as 2 —)■ 0, we choose gij = gij (which implies that fy = Tij). 
in these Fefferman-Graham coordinates then we have 


(E3) 

If we write both g^^, and 


5^1/ =9^i^ + z {Tij{x) - Tij{x)) dx^dx^ + O(z^). 


Plugging this into (E2) we see that 


r = y(r-r) + o(z6) 

f = - ^v*(r - r) j o(z®). 


(E4) 


(E5) 


where T = h^^Tij and Vi is the covariant derivative for h. We see the gauge condition is satisfied at this order if Tij is 
traceless and conserved. In the absence of a bulk gauge field, this is the case. However, the addition of a bulk gauge 
field means that, while the GET stress tensor is still traceless, it is no longer conserved. The boundary value of this 
bulk gauge held, Hi, is a source term in the GET which gives rise to a charge current Jb The expectation value of 
both this and the stress tensor are dehned through the variation of the effective action 

= J ^ (T,,) 6g^^ + (J^) S A. (E6) 

Demanding that this variation vanishes for coordinate transformations leads to the conservation law 

(E7) 


This means that in our case 


r = o(z®) 

c = - F,,F) + 0(z®). 


(E8) 


In fact, we’ve used the near horizon solutions as our reference metrics, and these have Fij = 0 and J* = 0. We can use 
this asymptotic form of in these coordinates to hnd an expansion for the transformation from Fefferman-Graham 
coordinates to the coordinates that satisfy the Harmonic gauge condition. It’s found to be given by 


x^ + 


z^ log z 


F,^F+0{z^). 


(E9) 


The presence of a logarithm in this transformation inevitably translates into a logarithm term in the boundary 
expansion of our solutions at fifth order. This means that they cannot be more than C^. 
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